Unraveling critical dynamics: The formation and evolution of topological textures 
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We study the formation of topological textures in a nonequilibrium phase transition of an over- 
damped classical 0(3) model in 2 + 1 dimensions. The phase transition is triggered through an 
external, time-dependent effective mass, parameterized by quench timescale r. When measured 
near the end of the transition (($^) — 0.9) the average texture separation Lsep and the average 



scale as Ls 



and Li, 



^o.46±o.04^ significantly larger than the single 



texture width L 

power-law ^freeze ~ r"'^^ predicted from the Kibble-Zurek mechanism. We show that Kibble-Zurek 
scaling is recovered at very early times but that by the end of the transition Lsepir) and Lw(t) 
result instead from a competition between the length scale determined at freeze-out and the ordering 
dynamics of a textured system. We offer a simple proposal for the dynamics of these length scales: 
Lsep{t) = ,fi(r) -I- Li{t - tfreeze) and L^(t) = ^2{r) + L2{t - tfreeze) whcre ^1 ~ t" and ^2 ~ 
are determined by the freeze-out mechanism. L\{t) — {^i)'^t3 and L2{t) = t^ are dynamical length 
scales previously known from phase ordering dynamics, and tfreeze is the freeze-out time. We find 
that Lsep{t) and Lw{t) fit closely to the length scales observed at the end of the transition and yield 
a = 0.24 ± 0.02 and (3 — 0.22 ± 0.07, in good agreement with the Kibble-Zurek mechanism. In the 
context of phase ordering these results suggest that the multiple length scales characteristic of the 
late-time ordering of a textured system derive from the critical dynamics of a single nonequilibrium 
correlation length. In the context of defect formation these results imply that significant evolution of 
the defect network can occur before the end of the phase transition. Therefore a quantitative under- 
standing of the defect network at the end of the phase transition generally requires an understanding 
of both critical dynamics and the interactions among topological defects. 
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I. INTRODUCTION 

We live in a world of broken symmetry. From the crys- 
talline structure of ice, no longer invariant under the full 
rotational and translational symmetry of liquid water, to 
the vacuum state of the standard model, the universe is 
filled with the patterns of symmetry breaking. Under- 
standing the dynamical processes by which symmetry- 
breaking phase transitions occur is therefore of substan- 
tial physical importance. Topological defects, forged in 
nonlinear interactions as the system seeks its true vac- 
uum, offer a window into the complex, nonequilibrium 
aspects of the phase transition. In systems with a topo- 
logically nontrivial vacuum manifold, fields in different 
spatial regions fall into different ground states and sym- 
metry breaking may be accompanied by the emergence of 
a well-defined network of topological defects at the end of 
the transition 0| . Simple topological defects include do- 
main walls, strings, monopoles and topological textures 
and are themselves the subject of much research in both 
cosmology and the laboratory Q |3|. Systems contain- 
ing topological textures are of particular interest as they 
provide an arena in which to explore myriad dynamical 
issues ranging from nonequilibrium phase transitions to 
the coarsening dynamics of quantum field theory. 



A. Topological defects in cosmology and the 
laboratory 

Very little is known about the exact nature of the early 
universe before nucleosynthesis. However, many par- 
ticle phenomenologies predict a unification of different 
forces at high energies, Tqut ^ lO^^GeV. If such unifi- 
cation exits then a rich pattern of symmetry-breaking 
phase transitions developed as the universe expanded 
and cooled. Phase transitions also occurred at the lower 
electroweak and QCD energy scales, although important 
physical details such as the exact order of the transitions 
remain unknown. 

In the course of symmetry breaking, causality con- 
strains the correlation length to grow no larger than the 
size of the horizon. Regions of space smaller than a hori- 
zon volume will rest in different vacua. For systems in 
which the vacuum manifold M has nontrivial homotopy 
groups n„(M), topological defects may exist where dif- 
ferent vacua meet. Topological defects formed in the 
early universe imprint themselves upon the fluctuations 
of the cosmic microwave background radiation (CMBR) 
and have been suggested as the generators of large scale 
structure formation (for a review see Much work 

has been done to compare measurements of the multi- 
pole moments of the CMBR anisotropics and numerical 
simulations of defect models. While the apparent lack 
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of agreement between theory and observations is not im- 
meadiately promising, the noisy nature of current mea- 
surements and the errors involved in doing large-scale 
numerical simulations of cosmological defect evolution 
make it difficult to rule out defect models of structure 
formation with any certainty (see for example |^). 

The study of topological defects in the early universe 
has been fruitful not only in attempting to explain pre- 
cise cosmological data, but also in influencing the direc- 
tion of cosmological research. One of the early moti- 
vations for the paradigm of inflation was the predicted 
overabundance of magnetic monopolcs If monopoles 
are formed in a phase transition before or during an in- 
flationary phase, their number density is safely diluted 
by the exponentially expanding scale factor. However, 
it has also recently been observed that in the very ener- 
getic process of reheating, as coherent oscillations of the 
inflaton decay into other fields, a nonthermal spectrum of 
topological defects may be produced, thus again raising 
the monopole problem |^]. The details of the reheating 
process and the possible formation of a defect distribu- 
tion at the end of an inflationary phase are currently an 
area of active study 0. 

Nonequilibrium phase transitions and the formation of 
topological defects are observed in many laboratory sys- 
tems and the study of topological defects has helped forge 
links between the study of high energy processes in the 
early universe and condensed matter systems 1^ . In liq- 
uid crystals, the formation of a network of topological 
string defects has been observed in the rapid cooling of 
the system through its critical point, demonstrating not 
only the viability of the Kibble mechanism in produc- 
ing topological defects but also the scaling of the defect 
network |9| . Other condensed matter systems such as su- 
perfluid '^e, have complex symmetry-breaking patterns 
and offer many potential analogs to early universe pro- 
cesses ||io| . 

Topological textures form in d spatial dimensions when 
the vacuum manifold M is such that Ild{M) is nontriv- 
ial. Unlike other topological defects, textured fields never 
leave the vacuum manifold but instead are "knots" of 
gradient energy, given relative stability by their topolog- 
ical winding around the vacuum. Like cosmic strings, 
the energy density in cosmological textures scales with 
the expansion of the universe so that textures formed in 
early universe phase transitions do not come to domi- 
nate the matter density. For this reason, global textures 
have been studied extensively as candidates for genera- 
tors of large-scale structure formation In con- 
densed matter, topological textures are studied in a wide 
variety of (mostly two-dimensional) systems, from super- 
conductors and quantum hall ferromagnets ||l4] to 
superfluid ^He-A 



B. Nonequilibrium phase transitions and topological 
textures 

The network of topological defects produced during a 
nonequilibrium phase transition contains dynamical in- 
formation about the evolution of the phase transition. 
A rapid phase transition can create a topological defect 
density far above equilibrium values. The defect density 
decreases as the system equilibrates after the transition. 
However the topological stability of defects allows any de- 
fect overdensity to last for a substantial period of time. 
Current understanding of the dynamics of topological de- 
fect formation is summarized by the Kibble-Zurek mecha- 
nism which describes the emergence of a single nonequi- 
librium length scale £, freeze characteristic of the initial 
defect density ||^. The freeze-out scale £, freeze is derived 
by combining phenomenological ideas of near-equilibrium 
critical dynamics and equilibrium critical scaling and pre- 
dicts a power-law dependence of the number of defects 
with the quench rate of the phase transition. Both the- 
oretical [p7|-pT| and experimental |^,^ effort has 
been focused on exploring the connection between the 
length scale of the initial defect distribution and the 
quench rate of the phase transition. These researches 
study the formation of topological defects that contain 
false vacuum in their cores: kinks in one dimension and 
vortices in 2 and 3 dimensions. Although (as explained 
below) the dynamics of a textured system can be qualita- 
tively different, the details of the formation of topological 
textures in a nonequilibrium phase transition are mostly 
unexplored (see however |l^). 

The dynamics of textured systems provides an ideal 
study for nonequilibrium phase transitions. Unlike other 
topological defects textures do not have a fixed size and 
the evolution of a textured system is complicated both 
by the poorly understood interactions between textures 
and by the changing scale of the texture itself. The rich- 
ness of texture interactions involves multiple dynamical 
length scales. In 2 -I- 1 dimensions it is observed that the 
phase-ordering of a system containing topological tex- 
tures is characterized at late times by at least three length 
scales: the (average) texture size, texture-texture sepa- 
ration and texture-antitexture separation [^,^ . In dis- 
tinction, systems of topological defects such as domain 
walls, vortices or monopoles are effectively single-scale 
and are described at late-times only by the average de- 
fect separation. In a nonequilibrium phase transition the 
multiple independent length scales of textured systems 
are determined by the phase transition dynamics. How- 
ever, the arguments of the Kibble-Zurek mechanism are 
ambiguous when applied to textured systems. Which, if 
any, of the multiple length scales observed in a texture 
distribution should we equate with the "frozen" correla- 
tion length? 
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C. Organization 

Section II reviews the properties of topological textures 
in 2 + 1 dimensions and the scaling violations observed 
in textured systems when quenched from a disordered 
phase. In section III, we present the Kibble-Zurek mecha- 
nism and derive the single nonequilibrium freeze-out scale 
^freeze characteristic of the defect network at the time of 
formation. Section IV contains the model of a nonequi- 
librium phase transition, the numerical techniques used 
in the evolution and a description of the evolution of the 
phase transition given by the numerical solution. We 
also define the length scales used to characterize the tex- 
ture distribution. Section V contains the main results of 
this work: an explanation of the quench rate dependence 
of the average texture separation and the average tex- 
ture width observed near the end of the phase transition. 
We show that the Kibble-Zurek mechanism is recovered 
at early times but that by the end of the phase transi- 
tion, formation dynamics and phase ordering dynamics 
are intrinsically linked. In Section VI these results are 
used to argue that textured systems carry an imprint of 
the freeze-out correlation length to late times and sug- 
gest this offers the possibility of novel, late-time mea- 
surements of dynamical critical phenomena. Concluding 
remarks are made in Section VII. 



divided into distinct topological sectors by winding num- 
ber n [p6| . In two spatial dimensions these windings are 
topological textures, topologically stable and static solu- 
tions with finite energy. For example, the following field 
configuration, a texture of negative unit winding, wraps 
around the vacuum 2-sphere exactly once. 



$1 



4ar sin((/)) 
r3 -I- 4a2 ' 



Aar cos((/)) 

J.2 ^ 4q,2 



4a^ 



4a2' 



(5) 



At the origin, <E> points in the — $3 direction while for 
r — > 00 it points in the opposite direction, $3. In be- 
tween, at r = 2a, the field points radially outward like 
a hedgehog. The parameter a characterizes the size of 
the texture and is independent of the parameters of the 
action Eq. The NLSM also possesses an exactly 

conserved topological charge Q, which is expressed as 
the integral over a topological charge density p, 



P = 1^$ • (5.* A dy$) 
47r 



Q 



(fxp{x,y) 



(6) 
(7) 



For the single winding texture, Eq. (^, the topological 
charge density has a particularly simple form 



II. TOPOLOGICAL TEXTURES IN TWO 
SPATIAL DIMENSIONS 

The properties of systems containing topological tex- 
tures differ in many ways from those containing other 
topological defects. To illustrate these differences, con- 
sider an 0(3) invariant model with classical action. 



and resulting equations of motion 

v'*« + A$,($'-m2) = 0. 



(1) 



(2) 



For the classical dynamics that is of interest here the 
value of the quartic self-coupling A and the mass parame- 
ter rn^ may be scaled to unity with an appropriate scaling 
of the field and spacetime units. 



t t\/\ni, 



(3) 



The vacuum manifold is characterized by <!>•$ = 1, topo- 
logically a 2-sphere. When the field is constrained to the 
vacuum the action is that of a nonlinear sigma model 
(NLSM). For asymptotically uniform fields, a vacuum 
field configuration is a map 



real space 



o2 

field space ' 



(4) 



1 

P= 



TT (r^ + Aa^y 



(8) 



with total topological charge Q = —1. The energy 
E = ^ J (Px(d^i ■ d^i) of the single winding configu- 
ration Eq. (g) is simply E = Stt, independent of the 
size a of the texture. In accordance with Derricks theo- 
rem [ p7| , the energy of a texture configuration in higher 
dimensions scales as a positive power of its size, demon- 
strating their instability to collapse. In lower dimensions 
textures are unstable to expansion. The stability of a 
texture configuration in any dimension can be altered by 
adding higher derivative terms to action Eq. (|l|), as was 
done in the Skyrme model of nucleons . In contrast 
to other topological defects such as strings or monopoles, 
the texture core size a is variable and independent of the 
fundamental parameters of the theory. 

Although the energy of a single isolated topological 
texture is independent of its size, systems with multiple 
textures are not static. Textured systems order under 
the constraint of conserved topological charge Q. How- 
ever, the details of texture interactions are not well un- 
derstood. It is observed that textures and antitextures 
can annihilate with each other and that more isolated 
textures can decay by unwinding |^,^. The details of 
the unwinding process depend on the model under con- 
sideration. For the action Eq.(|^) textures can decay by 
pulling the field off the vacuum manifold and unwinding. 
This is analogous to isolated textures in higher dimen- 
sions which are unstable to shrinking and unwind when 



3 



their gradient energy is large enough to move them off 
the vacuum manifold. In a NLSM where the fields must 
remain on the vacuum manifold, textures can unwind 
through higher-derivative terms present in the lattice dis- 
cretization. 



A. Topological textures and the violation of 
dynamical scaling 

When a disordered system is quenched into the ordered 
phase (e.g. by the removal of thermal fluctuations) the 
approach to a new equilibrium and associated long-range 
order is through a sequence of nonequilibrium states. A 
simple example is the Ising ferromagnet in 3 spatial di- 
mension. Above the Curie temperature, spins are ran- 
domly oriented and there is no net magnetization. A 
snapshot immeadiately following the quench of this sys- 
tem to zero temperature shows a system containing many 
domains in which the spins are coherently aligned. The 
net magnetization remains zero because the domains are 
randomly oriented with respect to each other. However, 
the true T = ground state of this system consists of only 
one domain and no domain walls. The system approaches 
this equilibrium through the dynamical process of coars- 
ening during which small domains shrink and large do- 
mains grow. The process of coarsening has been studied 
extensively in both condensed matter systems and in the 
early universe (for reviews see p9| and ||3C|]). The scaling 
hypothesis arose from these studies and states that, at 
late times, in a quench from a disordered state, the do- 
main distribution is statistically characterized by a single 
dynamical scale L{t) which grows in time. For simple dis- 
sipative systems, the approach to (T = 0) equilibrium is 
governed by model-A time-dependent Landau-Ginzburg 
(TDGL) dynamics for a nonconserved order parameter 



5$(a;,i) 
dt 



5^ 



(9) 



It has been shown both theoretically and experimentally 
that the average domain size L{t) grows with time as 



L{t) = (try- 



(10) 



which is also the power-law predicted by a dimensional 
analysis of Eq. (^). Not all systems coarsen in such a 
simple way and there is currently no general theoretical 
framework in which to explain why some systems scale 
while others do not. In particular, systems containing 
topological textures in one and two spatial dimensions 
are strong exceptions to the scaling hypothesis ||2^,^,^ . 
In a 2 -|- 1 dimensional 0(3) model with dissipative dy- 
namics and an instantaneous quench from the disordered 
phase, late-time coarsening of topological textures is de- 
scribed by at least three different length scales p^ , p5[ . 
These three length scales are the average defect-defect 
separation L^ep, the average defect width L^, and the 



average texture-antitexture separation Ltat- The scal- 
ing hypothesis is violated since at late times these three 
length scales are observed to grow with different powers 
of time 



-'sep{t) — £,ot^ 

L^{t)^^lt^ 

Ltat{t) 



t2 



(11) 

(12) 
(13) 



where is the correlation length in the disordered phase. 
It is remarkable that length scales describing the late- 
time ordering of textures retain a dependence on the ini- 
tial correlation length. In distinction, systems that obey 
the scaling hypothesis erase any memory of their initial 
state and at late times retain only a dynamical scale (Eq. 
(|l0| ) for the case of the ferromagnet). The scaling of sys- 
tems that contain topological defects such as vortices im- 
plies that the vortex distribution at late times is indepen- 
dent of the dynamical details of the phase transition. For 
textured systems it is not possible to have such a clean 
separation between the dynamics of the phase transition 
and the late-time texture distribution. In the context of 
a nonequilibrium phase transition is determined not 
from initial conditions but from the Kibble-Zurek mech- 
anism, thus connecting early-time critical dynamics and 
late-time coarsening. In Section VII we discus how this 
connection may provide late-time probes of nonequilib- 
rium critical dynamics. 



III. THE KIBBLE-ZUREK MECHANISM OF 
DEFECT FORMATION 

Before analyzing the dynamics of topological texture 
formation it is useful to review the basic Kibble-Zurek 
mechanism. The first estimate of the initial defect 
density in a cosmological context was made by Kib- 
ble |l|. The basic ingredients of the Kibble mechanism 
are causality and the Ginzburg temperature, Tq. The 
Ginzburg temperature is defined as the temperature at 
which thermal fluctuations contain just enough energy 
for correlated regions of the field to overcome the poten- 
tial energy barrier between inequivalent vacua, 



(14) 



where AF is the difference in free energy density between 
the true and false vacua and ^ is the equilibrium corre- 
lation length. In the Kibble mechanism, the length scale 
characterizing the initial defect network is set by the equi- 
librium correlation length of the field, evaluated at Ta- 
in a recent series of experiments |23 , the Kibble mecha- 
nism was tested in the laboratory. The results, while con- 
firming the production of defects in a symmetry-breaking 
phase transition, indicate that ^{Tg) does not set the 
characteristic length scale of the initial defect distribu- 
tion. These experiments were suggested by Zurek, who 
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criticized Kibble's use of equilibrium arguments and the 
Ginzburg temperature, Tq 

Combining equilibrium and non equilibrium ingredi- 
ents, Zurek offers a "freeze-out" proposal to estimate the 
initial density of defects. In the Zurek proposal, above 
the critical temperature, the field starts off in thermal 
equilibrium with a heat bath. As the temperature of 
the bath is lowered adiabatically, the field remains in lo- 
cal thermal equilibrium with the heat bath. Near the 
phase transition, the equilibrium correlation length and 
the equilibrium relaxation time of the field grow without 
bound as 



Tr = Tol e 



(15) 
(16) 



where e characterizes the proximity to the critical tem- 
perature, 



T,-T 



(17) 



and fj, and v are critical exponents appropriate for the 
theory under consideration. The quench is assumed to 
occur linearly in time 



t 

e — — 

T 



(18) 



so that for t < 0, the temperature of the heat bath is 
above the critical temperature and the critical temper- 
ature is reached at i = 0. The divergence of the equi- 
librium relaxation time as the heat bath approaches the 
critical temperature is known as critical slowing down. 
Critical slowing down results from the finite speed of 
propagation of perturbations of the order parameter. As 
the correlation length diverges, small perturbations of the 
order parameter {e.g., lowering of the temperature) take 
longer to propagate over correlated regions and therefore 
it takes longer to reach equilibrium. As the critical tem- 
perature is approached from above there comes a time 
I I during the quench when the time remaining before 
the transition equals the equilibrium relaxation time 



Tr{t*). 



(19) 



Beyond this point the correlation length can no longer 
adjust fast enough to follow the changing temperature 
of the bath. At time t* the dynamics of the correlation 
length "freezes" . The correlation length remains frozen 
until a time t* after the critical temperature is reached. 
In Zurek's proposal the correlation length at the "freeze- 
out" time t* sets the characteristic length scale for the 
initial defect network. Solving for the value of the cor- 
relation length at t*, the frozen correlation length, and 
therefore the initial defect density, scale with the quench 
rate as 



The power-law scaling of the topological defect den- 
sity with the quench rate has been verified using 
1- dimensional, 2-dimensional and 3-dimensional sim- 
ulations of phenomenological time-dependent Landau- 



Ginzburg equations 
turn field theory 



r7|-p0[ and in nonequilibrium quan- 
However, experimental efforts to 
test the Kibble-Zurek prediction in condensed matter sys- 
tems are inconclusive, lacking reliable error estimates and 
a mechanism to vary the quench rate p2| . In particular, 
quenches in '^He and high-temperature superconductors 
do not seem to produce defect densities at the level ex- 
pected from the freeze-out mechanism. 



IV. CRITICAL DYNAMICS 

In order to study the dynamical formation of topolog- 
ical texture the equations of motion Eq. (^ are modi- 
fied with an externally controlled time-dependent effec- 
tive mass. 



'^eff(^) = tanh( 



tc-t 



(21) 



where r is the quench rate and tc defines the critical point 
where spinodal instabilities, characteristic of a second or- 
der phase transition, begin to grow. Both in the early 
universe and in the laboratory, texture formation occurs 
in an environment containing many other degrees of free- 
dom. The effect of the environment is assumed to be that 
of a stochastic (white) noise driving term ^(x, t) and sim- 
ple ohmic dissipation, obeying a fluctuation-dissipation 
relation. The field equations of motion are 

- V'*» - 77*^ + «>.(*' - ^e//) = ^(a?, t), (22) 

where the fluctuation-dissipation relation is 

{C{x, t)^{x', t')) = 2Tf]5^{x - x')5{t - t'). (23) 

The use of a time-dependent mass in the equations of mo- 
tion of this 2-1-1 dimensional model deserves comment. 
At t = the effective mass is positive and the system 
fluctuates around the false vacuum |$| = 0. When t 
reaches tc the system is destabilized and the fields begin 
to roll to the true vacuum |$| = 1. This behavior mimics 
that of a phase transition in 3 -I- 1 spacetime dimensions. 
It is well-known that there is no continuous symmetry 
breaking and accompanying long range order in spatial 
dimensions of two or less 32 1. Therefore the time depen- 



at*) 



T 1+f . 



(20) 



dence of the effective mass does not come from changing a 
thermodynamic variable, such as temperature. However, 
though the 2-1-1 dimensional action of Eq. (|l|) has no 
finite-temperature equilibrium phase transition there can 
still be an effective dynamical transition. When energy 
is dissipated, a high temperature initial state dominated 
by short wavelength thermal fluctuations will eventually 
evolve into a state dominated by long wavelength topo- 
logical textures. Textures become quasi-stable excita- 
tions when the effective temperature drops substantially 
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below the texture energy of Stt. The use of an exter- 
nal time-dependent effective mass allows for a controlled 
passage between the fluctuation-dominated and texture- 
dominated phases. Within this controlled setting the role 
of the quench rate in the formation of texture can be 
isolated and studied effectively. The lessons learned are 
applicable to more realistic phase transitions, both in the 
early universe and in the laboratory. 



A. Numerical parameters and techniques 

The dynamics of the non-linear stochastic system 
was solved numerically. Eq. (|2^) was discretized on 
a 2-dimensional lattice using a standard second-order 
leapfrog method for the time evolution and a 

second-order spatial discretization for the laplacian y^, 



n 



n+l 



[jm 



1 



X 



n„[j][fc] 



dt 



(1 + x) 



(^(*n[j + l]W-4<I)„[j-][fc] 

+<g„[j-l][A;] + $„[j][fc + l] 



$„+i [j] [k] = [j] [k] + dffl„+i [j] [k] 



(24) 



where n = ^ and x ~ vdt- The labels [j][k] identify 
the lattice point and n labels the time step. At each 
time step and each lattice point the noise is generated 
through a sum of AI randomly distributed numbers pq] 



mk] 



g[j][fc] 24Af77r 
M ^ dx^dt ' 



(25) 



where —0.5 < 6'[j][A:] < 0.5 are a set of N"^ random num- 
bers. In the limit when Af — s- cx) the distribution 
approaches that of a Gaussian with the variance required 
for the fluctuation-dissipation relation Eq. (^) . In prac- 
tice M = 24 was used and no change in the dynamics 
was observed for larger M. The timestep was dt = 0.02 
and the lattice spacing was fixed at c?x = 0.5 on a square 
lattice with N = 1000 sites per side. There are three 
physical parameters in the system, the quench rate r, 
the dissipation constant 77, and the initial temperature 
T. The quench rate was varied between r = 10 and 
T = 150. Faster quenches (smaller r) deviate from the 
observed power-law behavior and asymptote to the be- 
havior characteristic of an instantaneous quench as has 
been previously observed Longer quenches are in 

principal possible but take substantially more computer 
time. The dissipation constant was chosen at ry = 1.0 to 
produce relatively overdamped behavior. The evolution 
equations were supplemented with initial conditions gen- 
erated by solving the equilibrium dynamics of Eq. (122 



with a value of the effective mass fixed at one quench 
timescale from the critical point, 



mljfit < 0) = tanh(l.O). 



(26) 



Run times for the system to reach an order parameter 
of (f) — 0.95 on a DEC 500 Mhz workstation varied from 
hours for r = 10 to days for r — 150. 



B. Properties of the texture distribution 

Field configurations obtained from numerical simula- 
tions provide an embarrassment of riches. In order to 
explore fully the participation of topological textures in 
the dynamics of the phase transition it is necessary to 
cull from the abundance of information contained in the 
snapshots of the topological field configurations, a few 
well-chosen functions whose time evolution provides a 
clear indication of the dynamics. The late-time coars- 
ening of a textured system is described by at least three 
different length scales |^,^, the average defect-defect 
separation L^ep, the average defect width L^;, and the 
average texture-antitexture separation Ltat- In the for- 
mation dynamics studied here we were able to measure 
reliably only Lggp and L^. These length scales are de- 
fined through the the topological charge density and the 
topological defect two-point correlation function. Con- 
sider the two-point defect density correlation function 



C,ir,t)^{\p{x;)\\pix- + r)\)-Mx)\y 

where () denotes average over the lattice, 
normalized correlation function 



Gir,t) 



C\ir,t) 



(27) 
Define the 



(28) 



Ci(0,t) 

G{r, t) is a statistical measure of the "lumps" in the de- 
fect density hmction |p(2^)|. The average texture width 
is proportional to the half-height scale ri of G{r, t) 
defined by 

G(ri,t) = 1/2. (29) 
To measure the average texture separation consider 



Qtotai = / d^x\p\. 



(30) 



Qtotai counts the total number of topological defects in 
the system without distinguishing between textures and 
antitextures. On dimensional grounds, 

Qtotai ~ J^^lys^ (31) 

Sep 

where Lsep measures the average separation of the topo- 
logical defects and L^yg is the physical size of the lat- 
tice which did not vary between runs. These two length 
scales, L^jSnd Lsep^ provide a detailed characterization 
of the topological texture network at the end of the tran- 
sition. In principle, each scale offers a different window 
into the noncquilibrium phase transition. 
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C. The dynamics of symmetry breaking 

Plots of the time evolution of the order parameter 
(f) = (<I'^) and the total number of topological defects 
Qtotai — J cPx I p{x,t) I are shown in figures and 
respectively. The plots correspond to an evolution 
with quench parameter = 10. In the symmetric state 
{t < 10.0) both (j) and Qtotai are close to zero. When 
t > 10.0 the external effective mass is negative and the 
fields, kicked by the small stochastic force £,{x,t), begin 
to roll to the true vacuum 0=1. The order parame- 
ter (j) is a measure of the number of lattice points where 
the fields have fallen to the true vacuum. As (p increases 
both the gradient energy and the number of defects in- 
crease. However as — *■ 1 most of the lattice is near 
the ground state vacuum and the defect formation pro- 
cess is over. After formation, the defect density does not 
remain constant but decreases slowly as texture-texture 
and texture-antitexture interactions dominate. 



V. UNRAVELING CRITICAL DYNAMICS 

Figures (^) and (jj) show Log-Log plots of the average 
number of textures and the average texture width versus 
the quench rate. The top graph in each figure displays the 
average number of textures and the average texture width 
measured near the end of the transition, taken to be the 
point where ($'^) — 0.9. The data are well-described by 
the power-law scalings 



T ^ _0.38±0.01 



^0.47±0.01 



-'Sep 



(32) 



The best-fit exponents were derived using equal weight- 
ing for all lattice measurements. As we found no reliable 
means to estimate the error in individual lattice mea- 
surements, the errors in the scaling exponents arise from 
statistical errors in the fit. Certainly, power-law scal- 
ing is expected from the freeze-out picture. However, 
the Kibble- Zurek mechanism also provides a precise pre- 
diction for the value of the power-law exponent as de- 
rived from equilibrium critical exponents. The model 
under consideration is a mean field model and the only 
role of the small fluctuations is to seed texture forma- 
tion. Therefore it is expected that = ^, the mean field 
value of equilibrium correlation length exponent. This 
expectation was verified by equilibrium lattice simula- 
tions. With ?7 = 1.0 and r > 10 the dynamics is over- 
damped at freeze-out. Therefore fj, = 2v = 1 and the 
Kibble-Zurek mechanism predicts the scaling ^(r) ~ r^, 
in obvious disagreement with Eq. (^2|). 

The larger observed exponent means that for longer 
quenches topological textures are both less numerous and 
bigger than the Kibble-Zurek mechanism predicts. How- 
ever, slower quenches take longer for the phase transition 
to complete since as r increases, the order parameter rolls 



more slowly. In fact, observations of the dynamics show 
that the time io.9 (measured from the critical point) it 
takes for the order parameter to reach ($^) = 0.9 in- 
creases with r in a simple way to. 9 ~ 1.5r. With more 
time to evolve between freeze-out and the end of the tran- 
sition, slower quenches allow for more interactions among 
the texture distribution. This suggests that the appar- 
ent power-laws observed at the end of the transition are 
not indicative of a single scale but a combination of the 
Kibble-Zurek mechanism and the dynamical length scales 
of the evolving texture network. 

Evolution of the power law exponents of the texture 
length scales is clearly indicated in the two lower graphs 
of figures ^ and |[ The lowest graph in each figure dis- 
plays the average number of textures and the average 
texture width measured very early in the transition when 
{^^) =0.1. At these early times the effect of coarsening 
dynamics is insignificant and power-law scaling indicative 
of the Kibble-Zurek mechanism is evident, 
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„0.22±0.01 



^0.14±0.01 



(33) 



Although the power law exponent of is slightly 
smaller than predicted by the Kibble-Zurek mechanism, 
the difference is not likely to be significant. Textures are 
not well-formed at these early times and the method for 
measuring L.^ is sensitive to the details of the two-point 
function. In the middle graphs of each figure, the num- 
ber of textures and the texture width were measured at 
intermediate times. Both show the crossover from forma- 
tion dynamics to coarsening dynamics as longer quenches 
develop a steeper power-law indicative of texture interac- 
tions. A similar picture of the evolution of Lg^p is given 
by the total gradient energy ^ 

Figures |3| and ^ provide clear evidence of the evolu- 
tion of the power-laws characterizing the length scales of 
the texture distribution. It is also possible to use the 
knowledge of texture coarsening dynamics to provide a 
quantitative understanding of the power-laws observed at 
the end of the transition. As discussed in section HA, the 
phase ordering of a textured system results in a dynam- 
ical length scale characterizing the separation of defects, 

Li{t) = ^0^3. The dynamics of this length scale for a 
fast quench (which enters the coarsening regime earlier) 
is shown in Figure ^. In the instantaneous quenches of 
p5| ^0 was identified with the initial correlation length. 
In a dynamical quench it is more natural to identify 
with the freeze-out correlation length. Combining the 
Kibble-Zurek mechanism with phase ordering dynamics 
I suggest that the evolution of the scale characterizing 
the separation of topological textures is given by 

Lsepit) = £, freeze + {S,freeze)^{t — t freeze ) ^ , (34) 

where ^freeze ~ is the length scale determined by the 
Kibble-Zurek mechanism, t freeze ~ t'^"' is the time from 
"freeze-out" and t is the time from the critical point. 
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Although Eq. (jsj) is shghtly complicated, its form is 
tightly constrained and the only free parameter is a, the 
exponent of the freeze-out correlation length. When t = 
to. 9 ^ 1.5t, the expected number of defects Q ^ j^— is 

Sep 

matched against the number of defects determined from 
the numerical simulations. A plot of the data and the 
best-fit length scale is shown in figure ^. To obtain the 
fit one data point is used to normalize the scale and a 
value of a is determined by minimizing x^- By using, in 
turn, each data point as a normalization a spread in the 
value of a is obtained. The best-fit value of a is taken 
as the average over all normalizations and the error is 
the standard deviation. As determined from these fits, 
a — 0.24 ± 0.02 in excellent agreement with a = 0.25 
expected from the Kibble-Zurek mechanism. 

The dynamics of L^,, the length scale characterizing 
the size of topological textures, is more complicated. Fig- 
ure ^ provides an example of texture width dynamics for 
T — 80. In the simulations the texture width was ob- 
served to grow in time with a power law near the simple 
coarsening dynamics predicted for an overdamped model 
L ~ ^2. However, this growth rate is different from 
previously reported studies of the late-time dynamics of 
the texture width. Consistent with previous studies ||2^ ] 
it is possible that L^{t) reaches its asymptotic dynam- 
ics only at later times. It is also possible that Lw{t) is 
contaminated at early times by topologically trivial spin 
configurations whose decay then dominates the dynam- 
ics. The resolution of this issue requires a more accurate 
method of counting textures e.g. by looking explicitly at 
the winding around the vacuum manifold. Whatever the 
mechanism. Figure |9| shows that Lw can evolve signifi- 
cantly between the early and late periods of the phase 
transition. As a concrete equation, I suggest that the 
evolution is given by 

Lw{t) — £, freeze + ~ t freeze) ^ • (35) 

freeze ^ IS the length scale determined by the Kibble- 
Zurek mechanism, t freeze ~ t"^^ is the time from "freeze- 
out" and t is the time from the critical point. As be- 
fore, when t = 1.5t, the expected width of the textures 
is matched against the width determined from the nu- 
merical simulations. A plot of the data and the best-fit 
length scale is shown in figure 0. The fits were obtained 
in the same way as for Lsep- As determined from these 
fits, (3 = 0.22 ± 0.07 also in agreement with (3 = 0.25 
expected from the Kibble-Zurek mechanism. 

It is clear from Figs, (j^) and (|^) that the power-laws of 
equation ( ^2|) can be explained by incorporating texture 
dynamics. Eqns. ^ and ^ offer a conservative, simple 
guess to the proper combination of formation and coars- 
ening dynamics. A more complete analysis, with better 
techniques for identifying textures, would fit Lsep{t) and 
(t) to functions with both general formation and coars- 
ening exponents. Nonetheless it is remarkable that these 
simple equations are in good agreement with the Kibble- 
Zurek mechanism. Even if the exact form of Eqns. js^ ) 



and (^) changes, the observation remains that the length 
scales characterizing the topological defect distribution 
at the end of the phase transition are a combination of 
defect interaction and formation dynamics. 

It is not surprising that defect interactions can influ- 
ence the properties of the defect network before the end of 
the transition. Although the long-wavelength dynamics 
is very slow until after the freeze-out time t freeze ~ 
(in the overdamped case), the time to complete the phase 
transition is tp.g ~ t allowing plenty of time for even 
partially formed defects to influence each other. In pre- 
vious numerical simulations in 1 and 2 dimensions [r^-[T9|] 
this effect was not noticeable (in overdamped evolutions) 
because the interactions between defects were extremely 
weak. In a recent simulation of global vortex formation in 
3 dimensions pO|] one might expect an observable change 
in the scaling exponents since the defect interactions are 
much stronger. However, the transition was at relatively 
high temperature and the defects were at least partially 
screened by thermal fluctuations. 

VI. TEXTURE COARSENING AND THE 
KIBBLE-ZUREK MECHANISM 

The close agreement between Eqns. (^) and ( |35| ) and 
the length scales observed at the end of the transition 
provides evidence that the a priori independent length 
scales, Lsep and L^, are formed by the same nonequi- 
librium dynamics of the Kibblc-Zurck mechanism. That 
this should be so is not obvious. In fact, the Kibble-Zurek 
mechanism can be questioned even in models which at 
late times depend only upon a single scale. In |Q it was 
argued argued that the length scales defined by the field 
correlation length and vortex density, in principle, derive 
from different attributes of the power spectrum. Only 
under restricted circumstances does the defect separa- 
tion follow the field correlation length and scale with the 
quench rate. However, the quench considered here was 
deep into the spinodal region and the early-time power 
spectrum is dominated by a single peak with momen- 
tum k = kmax in the low-momentum modes. When the 
power spectrum is strongly peaked, the freeze-out corre- 
lation length freeze ^ -j^ — determines a unique scale at 
early times to which all other length scales are connected. 

A. Late-time dynamics and Experiments 

The strong violation of the scaling hypothesis observed 
in the coarsening of 2 -I- 1 dimensional textured systems 
suggests the possibility of novel, late-time measurements 
of early-time critical dynamics. As discussed in section 
(IIA), the coarsening of a textured system with an instan- 
taneous quench retains "memory" of a length scale i the 
correlation length in the disordered phase. In a nonequi- 
librium phase transition it is the freeze-out correlation 
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length, determined by the Kibble-Zurek mechanism, and 
not the correlation length of the initial state, that scars 
the texture distribution at late times. For overdamped 
coarsening dynamics at late times, 



ML 



= constant — £,freeze- 



(36) 



The freeze-out correlation length can be measured by a 
late-time measurement of the (average) texture width 
and texture separation. It is the particular (and not 
well-understood) dynamics of texture-texture interac- 
tions that produces strong scaling violations and unusual 
memory effects. In contrast, experimental probes of vor- 
tex formation in nonequilibrium phase transitions of ^He 
and *He are hampered by the details of the interacting 
vortex distribution. Because these systems approach a 
scaling solution at late times, vortex interactions work to 
erase any memory of the initial state. This is a partic- 
ular problem in the '^He experiments where the vortex 
network is experimentally observable only at later times 
and the number of initial defects must be inferred us- 
ing detailed assumptions about the decay of the vortex 
tangle. 

Textured systems in 2 -I- 1 dimensions can be created 
in the laboratory. If their dynamics are approximately 
described by Eq. ( p^ ) then late-time experimental mea- 
surements of the Kibble-Zurek mechanism are possible. 
Topological textures can also be formed during phase 
transitions in 3 + 1 dimensions in superfluid '^He and 
the early universe. In 3-f 1 dimensions without additional 
higher derivative terms in the action, topological textures 
are unstable. In this case the dynamics of coarsening is 
complicated by texture collapse. However, as long as tex- 
tures exist in the system scaling violations and associated 
dynamical memory are possible. 



VII. CONCLUSION 

In this work, we exploited the multiple length scales 
of a textured system to study a nonequilibrium phase 
transition of an overdamped classical 0(3) model in 2 -I- 1 
dimensions. At very early times, we identified power-law 
scaling characteristic of the Kibble-Zurek mechanism and 
a single freeze-out scale. This suggests that the multiple 
length scales characteristic of the late-time ordering of a 
textured system derive from the critical dynamics of a 
single nonequilibrium correlation length. 

When observed near the end of the phase transition, we 
found the scaling of the texture separation Lsep (t) and 
the texture width L^ir) results instead from a compe- 
tition between the length scale determined at freeze-out 
and the ordering dynamics of a textured system. We ex- 
pect that this observation is not restricted to systems of 
topological textures or to low dimensions but is relevant 
anytime defect interactions are significant. It is not sur- 
prising that defect interactions can significantly modify 



the defect distribution even before the end of the phase 
transition. Well before defects are fully formed they frus- 
trate the system from its true minimum energy ground 
state and interact with eachother. 

The Kibble-Zurek mechanism provides a useful con- 
nection between critical dynamics and the length scales 
of the topological defect distribution observed at the end 
of the phase transition. However, power-law scaling pro- 
vides a rough and rather opaque window into the compli- 
cated nonequilibrium processes of the phase transition. 
As demonstrated explicitly here, power-law scaling can 
also arise from completely different mechanisms such as 
defect interactions. To understand the characteristics of 
the topological defect distribution at the end of the tran- 
sition it is therefore important to look closely both at the 
nonequilibrium dynamics of the phase transition and the 
evolution of the defect network. 
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FIG. 1. Plot of the order parameter (f){t) = (i$|^) versus 
time t for quench parameter = 10.0. Although the system 
is overdamped near freeze-out, the small oscillations in (j) near 
the end of the transition indicate that oscillations around the 
true vacuum are not completely overdamped. 




FIG. 2. Plot of the number of textures Q(t) = ^ifx\p{x,t)\ 
versus time t for quench parameter r = 10.0. As with the or- 
der parameter, the small "bump" indicates oscillations around 
the true vacuum near the end of the transition 
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FIG. 3. Plot of the average number of textures deter- 
mined from the topological charge Q ~ J d'^x\p{x,t)\ ver- 
sus quench parameter r for different values of the or- 
der parameter = ($^). Symbols denote lattice mea- 
surements and the two solid lines denote best-fit power 
laws. Early in the phase transition </> = 0.1 the best-fit 
power law is log(Q) = (1.31 ± 0.02) - (0.43 ± O.Ol)log(r), 
in good agreement with the Kibble-Zurek mechanism. 
At late times <j} — 0.9, the best-fit power law is, 
log(Q) = (3.17±0.03) - (0.77±0.02) log(r). Without texture 
interactions, the early and late-time power laws would be the 
same. However in this interacting system, texture unwinding 
and texture-antitexture annihilation steepen the power law 
observed at the end of the phase transition. The crossover is 
apparent at intermediate times (j> = 0.5 as defect interactions 
modify the number of textures for longer quenches. 
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FIG. 4. Plot of the average texture width versus quench 
parameter r for different values of the order parameter 
— ($'^). Symbols denote lattice measurements and 
the two solid lines denote best-fit power laws. Early in 
the phase transition = 0.1 the best fit power law is 
log(L„) = (0.42±)(0.14 ± 0.01) log(T). At late times = 0.9, 
log(L„) = (0.11±0.02)-(0.47±0.01) log(T). At intermediate 
times = 0.4, L„ is growing at a faster rate for the longer 
quenches. The early-time power law exponent is smaller than 
expected from the Kibble-Zurek mechanism. However, at 
early times the form of defect-defect correlation function. 
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FIG. 5. Plot of the average number of textures deter- 
mined from the gradient energy G = J Sx(d<^i ■ d^i) 
versus quench parameter r for difi'erent values of the order 
parameter (f> = ($^). Symbols denote lattice measurements 
and the two solid lines denote best-fit power laws. Early 
in the phase transition = 0.1 the best fit power law is 
log(G) = (2.34 ± 0.01) - (0.46 ± O.Ol)log(r), in agreement 
with the Kibble-Zurek mechanism. At late times (j) — 0.9, 
log(G) = (3.65 ± 0.10) - (0.94 ± 0.06) log(T). As in figure | 
a crossover is apparent at intermediate times (j> = 0.4. The 
late-time power law is slightly steeper than that determined 
from the topological charge because the gradient energy also 
includes topologically trivial configurations which decay more 
rapidly than textures. 
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FIG. 6. Plot of the average number of textures 
Q ~ J d^x\p{x) versus quench parameter r. For each 
quench the number of textures was measured at the point 
when {$^) = 0.9. Filled triangles denote lattice mea- 
surements. The solid line denotes the theoretically ex- 
pected number of textures computed from the length 

scale Lsep{t) — ^freeze + {^freeze)'^{t " t freeze)'^- The 

freeze-out correlation length determined from this fit is 
^freeze ^ ^o.24±o.o2^ excellent agreement with the Kib- 
ble-Zurek mechanism. 
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FIG. 7. Plot of the average texture width (in lattice units) 
versus quench parameter r. For each quench the average 
width was measured at the point when ($^) = 0.9. Filled 
triangles denote lattice measurements. The solid line denotes 
the theoretically expected width computed from the length 
scale Lw{t) = ^freeze + {t — tfreeze)^ ■ The freeze-out correla- 
tion length determined from this fit is (^freeze ~ t-0-22±o o7^ 
agreement with the Kibble-Zurek mechanism. 



FIG. 8. Plot of the average number of textures 
Q{t) — J cfx\p{x, t)\ versus time for a very fast quench t — 1. 
The solid line denotes data from a numerical simulation. The 
dashed line is the best-fit power law of the late-time coars- 
ening dynamics, Q{t) ~ t~^''^'^ . The coarsening exponent 
a = —0.72 agrees favorably with previous studies which found 
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FIG. 9. Log-Log plot of the texture width versus time 
for quench parameter r = 80. Time is measured from 
the critical point which occurs when tc = 80. Solid trian- 
gles denote data from a numerical simulation. The straight 
line is the best-fit power law of the coarsening dynamics, 
Lw{t — tc) ~ {t — tc)~^'^^. The coarsening exponent a = —0.62 
is close to the simple scaling exponent a = — \. 
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